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C^ . We study a generalization of the Weingarten model reduced to a point, which becomes 

the large- A'' reduced U{N) gauge theory in a special limit. We find that the [/(!)'' symmetries 
are broken one by one, and subsequently restored as [/(!)'' -^ U{1)''^^ ^ ■ ■ ■ —> U(l) -^ 
^\ ■ 1 -^ [/(I)'' as we change the coupling constants. In this model, we can develop an eflicient 

algorithm, and we can clearly see the phase structure of the large- A'' reduced model. We 
thus conclude that this model would be useful for the study of the unitary model. 

m 
> 

If^ • §1. Introduction 

o 
o, 

O : S = -(3NJ2^r (44^^^,) + KiV^Tr (4^^ - l)' , (M) 

Qh! where An {fi = 1,2,- ■ ■ ,d) are complex N x N matrices. This model is a generalized 






In this paper, we study a model defined by the action 



version of the reduced Weingarten model. 

The original Weingarten model^^ was proposed as a nonperturbative descrip- 
tion of the Nambu-Goto string. This model is defined as follows. Consider the 
/\ ' d-dimensional square lattice Z and introduce a complex N x N matrix Ax^^ for 

C^ ' each link connecting the sites x and x -\- fi in such a way that Ax+fi-fj, = ^x,/^- Then 

the action of the Weingarten model is given by 

d 
Sw = -NlY.T. Tr (^:.',m^x+a,^4+p,m4,.) + iV 5] 5] Tr (4,m^.,m) (1-2) 

X ^zfzv X fj,= l 

The partition function is given by 

Zw = dmNexp{—Sw), (1'3) 
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where the measure drriN is defined by 

dmiv = n n (^d [Re{A,,^)^j] d [Im{A^^^)ij]\ . (1-4) 



X,fJ. I,] 



Let Cj (i = 1, • • • , n) be closed contours on the lattice. Multiplying A^ along Cj and 
taking the trace, we obtain Wilson loops w{Ci). The correlator of 'w{Ci), defined by 

WiCi, ... ,Cn) = ^ f diUN exp {-Sw) ^^(Ci) • • • ]^w{Cn), (1-5) 

is evaluated as 

l^(Ci,---,CJ=iV2-2« Y, exp(-a(s)log7"^-/i(s)logiV2) (1-6) 

for large-A^, where S{{Ci}) is the set of surfaces on the lattice whose boundary is 
Ci U • • • U C„, a{s) is the area of the surface s, and h{s) is the number of handles 
of s. If we regard log 7"^ and -^ as the string tension and the string coupling, 
respectively, Eq. (jl-6p can be interpreted as the sum of random surfaces weighted by 
the Nambu-Goto action. 

Next, let us consider the reduced Weingarten model, ^^ whose action is given by 

d 

Srw = -Nj J2 Tr (^M^.44) + iV 2; Tr [a^A^) . (1-7) 

This action is invariant under the U{1)'^ transformation 

A^ ^ e^'^A^. (1-8) 

If this symmetry is not broken spontaneously in the large- A^ limit, the correlators 
of the Wilson loops of this model are equal to those of the original Weingarten 
model (|l-2p . Because the reduced Weingarten model has only d matrices, numerical 
calculations are more tractable. This model was studied numerically^^ in the cases 
d = 2 and 3, and it was shown that the Weingarten model does not describe smooth 
surfaces, but branched polymers. ^^ 

One possibility to overcome this difficulty is to consider the action (|l-ip . This 
action is motivated by the following observation. In the case of the Hermitian matrix 
model, we can describe a type OB string by flipping the sign of the double-well 
potential.^' Therefore, we expect that also in the case of the Weingarten model, 
worldsheet supersymmetry is introduced by modifying the potential, and it may 
prevent a worldsheet from falling into a branched polymer. 

Although there are several possibilities for modification, we choose Eq. (|l-ip be- 
cause it has the following interesting properties. First, at /3 = 0, this model consists 
of a set of d copies of a complex one-matrix model with a double-well potential. 
Secondly, at k = 00, the matrices ^^ are constrained to be unitary, and this model 
becomes the reduced U{N) gauge theory.^^ This model is worth studying in its 



own right, because of its relation to U{N) gauge theory and matrix models of su- 
perstring theory/" ^^ Thirdly, when /3 and k are finite, this model allows a lattice 
string interpretation similar to that allowed by the original Weingarten model, be- 
cause the relation (jl-6p also holds in this model, as long as no surface s intersects 
itself. Note that this action is the large-A^ reduction of the "interpolating model" 
proposed in Ref.,^^ in which the second term on the r.h.s. of Eq. (jl-2p is replaced 

with kN E, EUi ^ (4.,M^x,M - l)'- 

This model has been solved analytically only in the special cases (3 = 0^' and 
K = oo, d = 2.^'')'^^) We investigated the general parameter region numerically by 
Monte-Carlo simulation. 

The organization of this paper is as follows. In ^ we present the numerical 
results for d = 2, 3 and 4. We find d phase transitions that correspond to the 
partial breakdowns of U{1)'^ symmetry. Their A^ dependences and the existence of 
the hystereses indicate that they are of first order. At k = cxd, these phase transitions 
smoothly approach the known phase transitions of the large- iV reduced U{N) gauge 
theory.^^' In ^2.41 we show that at finite k, the U{1)'^ symmetry is restored if we 
further increase /3. Because in our model, numerical calculations are more tractable 
1*^1 . it is useful for the study of the unitary theory. ^ is devoted to conclusions and 
discussion. 

§2. Numerical results 

2.1. Numerical results for d = 4: 

We begin with the case d = 4. The phase diagram here is given by Fig. [TJ There, 
the symbols represent the points where the fluctuation of the action, {{AS)'^) = 
(5^) — (S)'^, diverges. The symbols + represent a third-order phase transition. Such 
a phase structure can be seen clearly only when N is sufficiently large, N > 30. 

At /? = 0, (jl-ip describes d copies of the complex one-matrix model, whose action 
is given by 

5onc-matrix = AtiVTr (^tA - l) . (2-1) 

The free energy is d times as large as that of the complex one-matrix model. The 
second derivative of the free energy with respect to n, C = K^-^{logZ), is given by 
Ref. 9): 

_C ^J #('^' + ('^-3)V^(A^ + 6) + i) ,(k<2) ^2.2) 

. (k > 2) 



iV2 d 

K 2 



The derivative of C with respect to k is finite but not continuous at k = 2, indicating 
a Gross- Witten-type phase transition. In Fig. [2] we plot C as a function of k~^ for 

*' We have applied the standard metropohs algorithm to each element of the N x N matrices. 
We have stored the matrices A^j^A^ in addition to the original A^'s. In this way, we can reduce the 
computational time for evaluating the variation of the action to order A'^, while it is proportional to 
N'^ in the case of the unitary matrix model. For details, see 3). 
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Fig. 1. Phase diagram of the d — 4 model. Phase transitions at points indicated by the symbols o 
seem to be of first order, while those indicated by + are of third order. Outside the solid line, 
representing /3k~^ — |, the action is not bounded from below. 



/3 = 0. Our numerical result accurately reproduces the analytic solution. We found 
a similar third-order phase transition at k~^ ~ 0.5 with < /3 < 0.15. For (3 > 0.15, 
we cannot see this phase transition clearly, because it is buried in the tail of the peak 
of the first-order phase transition. 

If we fix K at a sufficiently large value and vary /?, we find four critical points. 
We call them 01,(32,(3^ and P^ in ascending order. They converge to finite values as 
K — > oo. For K^^ > 1.0, these phase transitions seem to merge. 

Our result is consistent with that found in Refs. 12) and 13) at k = cxd. In 
Ref. 12) this parameter region was studied up to A^ = 16, and it was found that the 
order parameter P = ^ "^ iTr^d^p of the U{1)^ symmetry behaves asymptotically 

as P ~ 4 4= at large /?, where C ~ 2.162 is a fitting parameter. We plot P 

in Fig. [3l It is seen that at A^ = 16, our result agrees well with that of Ref. 12). 
Although the phase transition at Pi can be seen clearly, the phase transitions at 
/32,/33 and /?4 cannot be seen at this stage. For A^ > 30, we can see them clearly. 
This is consistent with the results of Ref. 13). 

In Fig. m we plot |TrA^| simultaneously for /z = 1, 2, 3, 4 at k = oo and A^ = 50. 
For /? < /?! ~ 0.190, we have |TrA^| = for any /i. This implies that C/(l)^ is 
not broken in this region. For Pi < j3 < P2 — 0.219, the quantity |TrA^| becomes 
nonzero. This implies that U{1)^ is broken to f7(l)^. In the same way, for P2 < P < 
Ps ~ 0.257, |TrA^| is nonzero for one of the A^, and hence U{1)'^ is broken to C/(l)^, 
for /^s < /3 < /34 ~ 0.307 |Tr^^| is nonzero for three A^ and hence C/(l)^ is broken to 
C/(l), and for P > P^, \TrA^\ is nonzero for all A^ and U{1)^ is broken completely. 
Similarly, for finite k, C/(1)^ breaks to ^7(1)^ at Pi, to U{lf at P2, to U{1) at /?3, 




Fig. 2. The "specific heat" C of the d — 4 modei at (3 = 0,N = 35. The solid curve represents the 
analytic result. 
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Fig. 3. Plot of P = -^ 5] iTrA^p for N = 16 a.nd N = 50. C ~ 2.162 is a fitting parameter. 
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versus /? for k = oo. Because 



and is broken completely at fS^. 

In Fig. [5] and Fig. [6l we plot P = ^ X^ 
there are hystereses around /?i and /32, these phase transitions are of first order. We 
conjecture that the other transitions are also of first order, although we have not 
seen yet clear hysteresis. The phase transitions at finite k also seem to be of first 
order. 

2.2. Numerical result for d = 3 

In this subsection, we consider the case d = 3. As we can see in Fig. [71 at 
large, fixed k there are three curves of first-order phase transitions. We call them 
Pi,/32 and /^s in ascending order. U{1)^ is broken to ?7(1)^ at /3i, to U{1) at P2, and 
completely at P^. The quantities /3i,/32 and (3^ converge to finite values as k"^ — > 0. 
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Fig. 4. (7(1)* breakdown in the unitary matrix model at A'^ = 50. The quantities i|Trj4^| {fi = 
1, 2, 3, 4) are plotted simultaneously. 
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Fig. 5. Hysteresis in the phase transition of the 
unitary matrix model around /3i. N = 50. 



Fig. 6. Hysteresis in the phase transition of the 
unitary matrix model around (32- N = 50. 
^=i7.EJTrA,r 



In Fig. [8] we plot |Tr^^| simultaneously for /i = 1, 2, 3 at k,~^ = 0. For /3 < /3i ~ 
0.30, |Tr^^| = for all /x. This implies that U{Vf is not broken in this region. For 
/3i < j3 < 132 — 0.40, |Tr^^| becomes nonzero for one A^^. This implies that U{1)^ 
is broken to U{1)'^ . In the same way, for /32 < /3 < /Ss ~ 0.57, |Tr^^| is nonzero for 
two A^, and hence C/(l)^ is broken to C/(l), and for /3 > Ps |Tr^^| is nonzero for all 
A^, and C7(l)^ is broken completely. For k~^ > 1.5, the three curves of the phase 
transitions seem to merge. 



2.3. Numerical results of d = 2 

In this subsection we consider the case d 



2. In this case, if k 



= 0, this model 

2 92 



is equivalent to the Gross- Witten model, and the specific heat C = (3 -B^Z is given 
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Fig. 7. Phase diagram of the d = 3 model. Outside the sohd curve, the action /9k ^ 
bounded from below. 
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Fig. 8. U{iy^ breakdown in the unitary matrix model at Af = 50. -^iTrA^jl (/i — 1, 2, 3) are plotted 
simultaneously. 
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A third-order phase transition takes place at /3 = ^. This can also be observed by 
numerical analysis (see Fig. [9|). 
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The phase diagram for A^ = 50 is given in Fig. 1101*^1 At large, fixed k there 
are two curves of first-order phase transitions. We call them Pi and P2 in ascending 
order. They correspond to the breakdown of U{1)'^ symmetry. If we increase P with 
K fixed, first U{1)'^ is broken to ^(1) at Pi, and it is broken completely at 02- 

At small K, Pi and P2 seem to merge. They seem to go to oo as k~^ — > 0. This 
is consistent with the analytic result, in which C/(l)^ is not broken. 

In Fig.[Tn]there is a line of the symbols o. Beyond this line, the U{1)'^ symmetry 
is restored. We discuss in detail this novel phenomena in the next subsection. 
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Fig. 9. The specific lieat of tlie d = 2 model at k~^ = 0,N = 25. The solid curve is the analytic 
result (EH). 



2.4. Restoration ofU{l) symmetry 

In this subsection we show that the U{1)'^ symmetry is restored as k~^ or 
P becomes larger. This restoration takes place in the well-defined region where the 
action is bounded from below. 

In order to see whether the action is bounded from below, it is enough to consider 
the quartic term 



d d 

Iquartic = -pN ^ Tr (^^t ^M^-) + «A^ J] Tr (^A^A^AlA, 



(2-4) 



This is positive definite if and only [i[^ Pk ^ < ^^. In Fig. \TT\ we plot P 



'W S/i I^t^mI^ at k ^ = 2.0. For P ~ 0.315 and 0.340, there are discontinuities in 
the value of P. These correspond to the breakdowns of U{1)'^ symmetry. As we 



*■* For K — 2 and 0.25 < /3 < 0.4, we studied up to A*' = 75 and found no other phase transition. 
**' This condition comes from the following steps. First, the inequality of the arithmetic and 
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Fig. 10. Phase diagram of the d — 2 model. Outside the solid curve, representing /3k~^ = 1, the 
action is not bounded from below. 



f+++++++++++++++++++++++++++++++++++++++ 



approach the boundary of the well-defined region, f3 = 0.5, the value of P decreases 
to zero, which suggests that the U{1)'^ symmetry is restored. In FigllOl we denote 
this phase transition by the symbols o. 

Although the signal is not clear, we conjecture that the restoration of f/(l) 
takes place also for d = 3 and 4. Because such a phase transition would take place 
near the boundary of the well-defined region, we believe that it does not occur in 
the unitary model. 

§3. Conclusions and discussion 

We have found d first-order phase transitions corresponding to one-by-one break- 
downs of the [/(l)"^ symmetry. As k~^ — > 0, they are smoothly connected to those 



geometric means, 



implies 



2Re Tr AB^ < Tr AA^ + Tr BB 



2Re Tr (^A^A^aIaIJ < Tr (^AIA^A^AIj + (^^^yljl^U- 
<■ ir \A^Ai_iA^Ai_i j + [A^jAi^Aj^Ai^j 
Thus, summing the spacetime subscripts, we get 

d ^2 

J2 Tr [aIaIA^A^) < (rf - 1) ^ Tr (^A^) . 
Equality here is reahzed, e.g., in the case of the unit matrix. 
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Fig. 11. Plot oi P= ^Y.^ iTrA^I^ for k"^ = 2.0 and A^ = 50. U{lf is broken at /3 ~ 0.315 and 
13 ~ 0.340, and then restored at /? ~ 0.48. 



found in Ref. 13). We also found a first-order phase transition corresponding to the 
restoration of U{1)'^. These phase transitions can be clearly seen for A^ > 30. In 
the region where c (< d) of the U{1) are not broken, our model can be regarded 
as a model defined on a c-dimensional lattice. In particularly, at k~^ = 0, it may 
represent a c-dimensional Yang- Mills theory coupled to {d — c) matter fields. 

Although the first-order phase transitions mentioned above may not allow a 
continuum limit, there is a continuous curve of third-order phase transitions which 
may be equivalent to a lattice string theory. In order to determine whether it has 
a continuum limit, we must study the expectation values of Wilson loops and the 
Creutz ratio. It would also be interesting to study the "interpolating model" ,^' which 
does not reduce to a point. Because our reduced model (jl-ip and the interpolating 
model are not equivalent for (5 > f3i, where C/(l) is broken, there is still a possibility 
that the interpolating model describes a continuum string in this region. 

In the large- A^ U{N) gauge theory on a lattice, there is a bulk phase transition, 
which does not involve the breakdown of any symmetry. In our simulation, however, 
we do not find a corresponding phase transition. However, this is not a contradiction, 
because in the large-A^ reduced U{N) gauge theory, the bulk phase transition is 
hidden by the breakdown of C/(l) symmetry. ^'^^ We hypothesize that a similar 
phenomenon occurs also for k < cx). 

The large- A^ reduced U{N) gauge theory can be regarded as a toroidal compact- 
ification of the bosonic part of IIB matrix model. ^^ In this interpretation, the C/(l) 
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syinmetry bocomes the translational symmetry, and there emerge spacetimes with 
various dimensions, depending on the coupling constant. 
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